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Abstract 

We present a simple stochastic algorithm for generating multi- 
plicative processes with multiscaling both in space and in time. With 
this algorithm we are able to reproduce a synthetic signal with the 
same space and time correlation as the one coming from shell models 
for turbulence and the one coming from a turbulent velocity field in 
a quasi-Lagrangian reference frame. 
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1 Introduction 

The multifractal language for turbulent flows has been introduced about 
20 years ago in order to describe the anomalous scaling properties of turbu- 
lence at large Reynolds numbers jl], ^| . Beside any particular interpretation, 
the multifractal formalism exploit the scale invariance of the Navier-Stokes 
equation by taking into account fluctuations of the scaling exponents. To 
be more quantitative let us consider the Navier-Stokes equations: 

dtu + u-Vu— — Vp + i/Au (1) 
P 

where u is the velocity field describing a (homogeneous and isotropic) tur- 
bulent flow. For v = the Navier-Stokes equations are invariant with 
respect to the scale transformation: 

r -> Xr \ h u t -> Xt 1 -' 1 (2) 

Then, following Kolmogorov, it is assumed that at large Reynolds numbers 
[y — > 0) the rate of energy dissipation is constant. As a consequence, 
h = 1/3, if no fluctuation on h are present. The above reinterpreta- 
tion of the Kolmogorov theory naturally opens the way to describe in- 
termittent fluctuations in turbulent flows. Following the original idea by 
Parisi and Frisch, many possible values of h are allowed in turbulent flows. 
Each fluctuation h at scale r is weighted with a probability distribution 
p h (r) ~ r 3-DW 

Since its first formulation, the multifractal model of intermittency have 
been applied to explain many statistical features of intermittency in a uni- 
fied approach. The final goal of many theoretical investigation is to com- 
pute the function D(h) starting by the equation of motions. In some simple 
although highly non trivial case such a goal has been recently reached for 
the case of the Kraichnan model of a passive scalar || . 
One of the key issue in the multifractal language of turbulence is to un- 
derstand in a more constructive way what is a multifractal field and how 
the fluctuations of h are related to the dynamics of the system. In order 
to develop any systematic theory for computing D(h) starting from the 
equation of motions, one has to handle a complex non linear problem: the 
way in which a perturbative scheme may be developed strongly depend on 
a reasonable ansatz on the time-space properties of the probability distri- 
bution. It is therefore crucial to understand how we can formulate the most 
general form of multifractal random field which is consistent with the time 
and space scaling properties of the Navier-Stokes equations. 
One possible interpretation of the multifractal formalism is to observe that 
for any r < R, the multifractal theory predicts: 



Su(r) = W(r, R)Su(R) 5u(r) = u(x + r) - u[x) (3) 
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Then, according to the scaling properties of u, the quantity W(r, R) is a 
random quantity proportional to {-^) h ■ It turns out that for r\ < r 2 < r$ 
we have 



Equation (|J) tells us that one possible interpretation of multifractal field 
is to assume that fluctuations at scale r are described by a random mul- 
tiplicative process. The random multiplicative process is also somehow a 
simple way to mimic the energy cascade in turbulence. Actually, a general 
formulation of multifractal random fields based on random multiplicative 
process was first presented in ||, |j| by using a wavelet decomposition of 



One obvious limitation of random multiplicative process is the absence of 
any time dynamics in the field, as one can immediately highlight by consid- 
ering space-time correlations. Space-time scaling is a crucial and delicate 
issue when considering multifractal fields for the Navier-Stokes equations 
|U §]• It is the aim of this paper to understand how one can exhibits a 
multifractal field whose space and time scaling is consistent with the scal- 
ing constraints imposed by the Navier-Stokes equations. In Section 2 we 
introduce the technical problem shortly reviewed in this introduction by 
using a rather simplified language. In Section 3 we discuss several implica- 
tions of the results obtained in Section 2, with a particular emphasis on the 
consequences for the fusion rules as introduced in |^|, 0. In Section 4 we 
outline our conclusions and we discuss future extensions of our research. 

2 Multi-scale and multi-time stochastic sig- 
nals 

To simplify even further our argument we can concentrate an a typical 
fluctuation at a given scale, i.e. disregarding space position. We introduce 
the scale hierarchy l n — Iq ■ A~ n , in terms of the scale separation A > 1, and 
the velocity differences w n — v(x + l n ) — v(x). We assume that the scaling 
properties of w n are consistent with the dimensional constraints imposes 
by the Navier-Stokes equations, i.e. 



W(r u r 3 ) = W(r u r 2 )-W(r 2 ,r 3 ) 



(4) 



the field. 



d t w, 




(5) 



If w n shows multifractal scaling, we may write: 





where Q(p) is a non linear function of p and (. . .) is an average over the large 
scale statistics. Following Parisi and Frisch, we know that the multifractal 
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scaling (^) can be derived by assuming that w n ~ 1% with probability 
ll- D(h \i.e. 

«)~ fdhlt + "- D{h) (7) 



Indeed by means of a saddle point evaluation of the previous integral, one 
obtain the explicit expression for Q p in terms of D(h): 

C(p) = inf [ph + 3 - D{h)] (8) 

h 

Supposing one wants to keep into account also the time correlations, the 
constraint (|^) implies that 

C m (t) = (w n (tyw n (0)i) ~ fdhl^+^+ 3 - D ^f p , q (T/r n ) r n = l n /w n 

(9) 

where r n is a random time (the eddy turnover time) and the functions f p , q 
are dictated by the dynamic equations. Expression (||) has been introduced 
in H and analyzed in details in ||. We underline that, as a consequence of 
(||), we can predict the scaling properties of quantities like d ^'^} T \ We 
now want to understand how to define a random process satisfying both 
multiscaling in space (^|) and multiscaling in time (^|). In a more general 
way, we would like to exhibit random multifractal fields with prescribed 
dynamical scaling. It is known that the multifractal scaling (0) can be ob- 
served for random multiplicative process. Let us introduce the (positive) 
random variable A n and let us indicate with P(A n ) the probability distri- 
bution of A n . Then, by defining w n — (11™= l A*) w o an( ^ by assuming that 
the random variables Ai are independent, one obtains: 

(n \ P n 
llPiAJdAi = (AT = Co(p ) = l og ((^))/logA 

i=l / i=l 

. (10) 
We want to generalize expression (|10|) in order to satisfy the dynamical 

constrain (||). At each scale l n we introduce the random time r„ = l n /w n . 

The generation of our signal proceeds as follows, we extract A n with 

probability P(A n ) and keep it constant for a time interval [t, t + T n ]. Thus, 

for each scale l n , we introduce a time dependent random process A n {t) 

which is piece- wise constant for a random time intervals [<4 , tn + T n ] , if 

(k) 

tn is the time of the fcth jump at scale n. The corresponding velocity field 
at scale n, in the time interval tn < t < t n k ^ + r„, is given by the simple 
multiplicative rule: 

w n {t) =i4 n (t)w n _i(4* ) ) 

What is important to notice is that at each jumping time, t„\ t„ \ ..,t!£\. 
only the local velocity field is updated, i.e. information across different 
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scales propagates with a finite speed. This is one possible and relatively 
simple way to take into account the constraint ([5]) . 

To give a visual idea of how the algorithm works we show in Figure [l] 
the time behaviour of the regeneration times for several scales. It is evident 
that there are short time-lag where the chain of multipliers is not given by 
an exact multiplicative process (this happens every time a small scale has 
to be regenerated but the ancestors are not yet dead). 

In Figure || we show the scaling behaviour of the third order structure 
functions Ss(ln) = (wy obtained by a numerical simulation of a time de- 
pendent random multiplicative process with P(A) = p a 8(A — A a )+pbS(A — 
Ab), where A a = 0.2, Ab — 0.6, pb = 1 — p a and p a has been chosen such 
that £(3) = 1. Although S3 shows a very well defined scaling, the value 
of C(3) is greater than what is predicted by ( |lC| ) (in Figure || the slope — 1 
is shown for comparison). This effect shows that the "real space" scaling 
(iv?) is renormalized by the presence of the non trivial time dynamics of 
the multipliers. 

In order to understand why ( |To| ) cannot be used to predict the scaling 
exponents, let us understand which is the effect of the time dynamics for 
a given scale l n by assuming that at scale l n -\ the variable w n -\ is kept 
constant. Let T be the time used for time-averaged quantities and let N a 
and Nb be the number of events where the random variable A is equal to A a 
and Ab respectively. We next introduce the quantities r a = l n /(A a w n -i) 
and Tb — l n /(AbW n -i), the times associated to A a and Ab. By using 
our definition we can write: N a r a + NbTb = T, N a = p a N, Nb = 
p b N, N a + N b = N. It then follows: 

«> = K-r) \ J dtA(ty = {wl_ x ) ±{r a N a Al + r b N b Al) (11) 

The above expression can be further simplified and we finally obtain: 

™-<<*> T&££ (12) 



The consequence of (12) is that the scaling exponents C(p) are renormalized 



according to the expression: 

C*(p) = Co(p-l)-Co(-l) (13) 

where the number Co(p) are the "bare" scaling exponents, i.e. those com- 
puted by using P(A) according to (|l(i|). 

Expression (|l3|) have been obtained by using the simplified assumption 
w n -i = const. In the general case, i.e. all variables w n are fluctuating, one 
needs to generalize the above discussion. A possible way is to write: 

Cfl(p) = Co(p-«)-Co(-a) (14) 
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where the number a (not necessarily integer) depends on the details of 
P(A). 

We have checked expression (|l4|) for a number of different choices of 
P(A). Here we present the results for P(A) being a log-normal distribution, 
i.e. for 

Co(p)=p/io-^V (15) 

By using ( |l4| ) we obtain: 

C R {p)^p{K + c t a 2 )-^a 2 p 2 (16) 

In Figure (^) we show Cr(3) obtained by a set of numerical simulations for 
different values of a and ho chosen in such a way that Co (3) = 1- In this 
case ( |l6| ) can be written as Cr(3) = 1 + 3atr 2 . As one can easily observe, 
the prediction of (^|) is verified with very good accuracy with a value of 
a close to 2. In Figure (^) we show the value of Cr(p) for p = 2, ..,6 as 
obtained by direct numerical simulation for a = 0.03. The dashed line 
represent the estimate ( |l6| ) with a = 2: a very good agreement is observed. 

The above discussion can be generalized for multifractal fields and for 
any particular choice of P{A). In all cases, a re- normalization of the scaling 
exponents, as predicted by ([l6|), should be expected. At the same time, a 
non trivial time correlation is introduced for the variables w n . In Figure 
(||) we show (w n (T)w n (0)) and (w n +6(T)w n {0)) for n = 3, obtained by a 
numerical simulation of the time dependent random multiplicative process 
with a log-normal distribution with parameters a = 0.03, Cr(3) = 1- As 
expected, the correlation (w n +e(T)w n (0)) increases for small r and then 
goes to 0. The pick at a time lag larger than zero is due to the presence of 
a non trivial time dynamics, i.e. multipliers at different scales need some 
time to realize that their ancestor have changed their status. This is meant 
to mimic the non-trivial time dynamics of the turbulent energy transfer. 
This behaviour is also observed in the numerical simulation of deterministic 
shell models as reported in j|. 

Finally let us check whether the quantities w n satisfies the scaling con- 
strain imposed by (||). We first observe that the correlation function 
(w n (t + T)w n (t)) goes as exp(— _B|r|). This is due to the fact w n as a 
function of time is not differentiable. In order to check whether (||) is sat- 
isfied, we observe that B ~ k n w^ ~ const if Cr(3) = 1. In Figure (||) we 
plot the timescale B as a function of n obtained by a time dependent ran- 
dom multiplicative process, using a log normal distribution with a — 0.03, 
Co (3) = 0.83 and Cr(3) = 1- As one can see, B is fairly constant in the 
inertial range. 
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3 Numerical results 

The re-normalization effects discussed in the previous Section can be fur- 
ther investigated by considering the case of a passive scalar. In this case 
the constraint (||) should be written as: 

d9 n ~0j^ (17) 

^n 

where 9 n = 9(x+l n ) — 9(x) in analogy with the definition of w n . We assume 
that (6P n ) ~ ln (p) . Let us assume that a suitable representation of 9 n is 
given by a time dependent random multiplicative process as described in 
Section |^. In this case, however, the random time r„ is not correlated with 
the value of 9 n or the related random multiplicative variables. Therefore, 
we should not expect any re-normalization for the scaling exponents H{p). 
This is indeed the case as shown in Figure (Q), where we plot the scaling 
of (9 n ') for a log-normal random multiplicative process with <r = 0.03 and 
ho such that H(3) = 1. The updating times have been chosen with an 
independent random distribution mimicking the velocity fluctuations, r„ ~ 
l n /w n The dashed line corresponds to a slope —1. 

Our definition of time dependent random multiplicative process can be 
very useful in investigating the behaviour of the so called fusion rules for 
the quantities like (w^ +m w^). Following [|[ 0. we can write: 

« +m <) = C p , q (m)^h^- « +q ) (18) 

where C Pig is a constant for large m. Actually a direct measurements of 
(w n+m w ( l l ^ in turbulent flows at high Reynolds numbers and in direct nu- 
merical simulations, confirm the validity of (|l^) with C VA < 1 for large m 
[ p"3| . It is interesting to observe that our time dependent random multi- 
plicative process satisfies fusion rules with C p . q < 1. It is indeed possible 
to show this result by the following argument. Let us consider two scales 
/„ and l n +m- For fixed time, the quantity w n and w n + m are not necessary 
product of the same random variables. They will feel the same chain of 
multipliers only for some larger scale l n > with n' < n, while for smaller 
scales both chain of multipliers leading to the instantaneous value of w n 
and w n + m may have changed. Thus we should have: 

In' \ I n \ ln+m \ 

«< +m ) = ( n a i +p ){ii A i)(ii A i) 

U-l / \k=n> I \k=n' I 

where is the random multiplier acting between scale k — 1 and scale k, 
i.e. Wk+i = AkWk- The above expression gives the following result for the 
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compensated fusion rules: 



< +m )(< +q ) (AP) n+m (AP+i) n m 



that is: 



(« +m ) «) _ ( (Ai+n 



< 1 (21) 



Let us notice that for n' = n, which corresponds to random multiplicative 
process without time dynamics, the r.h.s of the above expression is just 1. 
The above equation should be considered as a qualitative prediction. In 
general we expect n' to be a function of both p and q. 

We have checked our prediction by several simulations for different 
choices of P(A). In Figure (||) we show the quantity Ci.2(jn) defined as in 
( p"8| ) for n = 3 and for m = — 1, ..11. As expected Ci.2(m) = 1 for m = 
and for m > is a slowing decaying function of m which reaches a plateau 
only for very large m. 

As a consequence of our analysis, we can also predict that the quantity 
C p ^ q (t), defined through: 

( w p \ 

(w n+m (t)p Wn (oy) = c p , q (t y /"+ n/ m » i 

should increase with t. This is indeed the case as one can observe in Figure 

(I)- 



According to (21), the asymptotic value of C p . q (m) for large m is a 
function of the intermittency, i.e. C Pj<2 (oo) becomes smaller for larger value 
of the anomaly (r.{p) + Cr{q) — (r(P + The qualitative prediction of 
has been checked against numerical simulations of time dependent 
random multiplicative process. In Figure ([l0|) we plot the asymptotic value 
of (^2,2(00) = ^"im^ooC^, 2(w) for a log normal P(A) as a function of 
Cr(4) — 2Cr(2). As one can see the qualitative prediction is confirmed. 

A further inspection of the numerical simulations reveals that C Pi9 can 
be written as 

C M (m) = l-A M /(m) (22) 

where f(0) — and /(oo) = 1. While A Piq is a function of (r(p) + 
Cr{q) ~ Cr(p + 9)1 we may wonder whether the function /(m) is somehow 
universal. Although a definite conclusion cannot be reached by looking at 
the numerical simulations, still our results seem to indicate that /(to) is 
either universal or is a function weakly dependent on intermittency. This 



can be seen from Figure (11) where we plot /(to) for different values of 
Ch(4) - 2Cr(2) . 
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We already observed at the beginning of this Section that a time de- 
pendent multiplicative process for a passive scalar does not show a re- 
normalization of the scaling exponents H(p). However, the argument used 
to derive ( f2l| ) can be applied even if the random times r„ are not correlated 
to the multiplicative process, as in the case of the passive scalar. In Figure 
T2|) , we show the quantity G(m) defined by the relation: 

(m +m ) = G(m)^^ (6t) (23) 

for the passive scalar and m = — 1..11. As expected, the fusion rules are 
satisfied only asymptotically with prefactors G(m) < 1. 

Let us also notice that the coefficients C Pi9 of the fusion rules are not 
fully defined, in terms of the time dependent random multiplicative process. 
Let us consider the variable 

U n = g n W n 

where g n is a random variable, independent of n, with the same probability 
distribution for any n, i.e. {g ni 9n 2 --9n k ) = (g ni ) (flO •• (flO and p (9n) 
does not depend on n. It follows that the fusion rules for u n satisfy: 



«+m<) = (9 P ) (9 q ) «+mO 



which gives 



(u p \ 1 

« +m <) = (g p ) (g q ) C P , g ^f « +? ) 



(g p+q ) 



Note that, for each n, w n is the time dependent multiplicative process and 
that the scaling properties of u n and w n are the same. The above equation 
implies that the fusion rules coefficients C Pjl? for u n becomes: 



(g p ) (g q ) 

(gP+1) 



C p> g -> C P}<] w n -> g n w n (24) 



Equation ( |24| ) implies that the asymptotic value of C P:q for large m is 
fixed up to a number (less than 1) linked to a scale invariant probability 
distribution function. 

Before closing this Section we would like to notice that the renormalized 
exponents Cr(p) are associated to a less intermittent field with respect to 
the bare exponents Co(p)- This can easily be deduced by using equation 
([Til). Thus we can predict that whenever the random times r„ are not 
correlated to w n we should observe (for a given P(A)) an increase of inter- 
mittency. One may wonder whether this qualitative prediction may have 
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any experimental evidence. This is indeed the case. Let us consider a shear 
flow. It has been noted that whenever the mean shear S is large enough 
intermittency increases. For large scale the characteristic time scale of the 
dynamics should be dominated by the shear effect and we would expect the 
characteristic time scale independent on the scale [|l0| [ll], [l2|]. In this case, 
therefore, if we describe intermittency as a random multiplicative process, 
time dynamics does not lead to a re-normalization of the scaling exponents, 
i.e. we should observe the scaling exponents £o(p) < Cr{p)- It is suggestive 
to think that the increase of intermittency in shear dominated flows may be 
understood in terms of the absence of re-normalization in time dependent 
multiplicative process. 

4 Conclusions and discussions 

We have introduced a simple multiplicative process which embeds intermit- 
tency both in time and in space. This allow us to generate a signal which 
respect the constraint (imposed by the Navier-Stokes equations) given by 
equation <^): dtiv n ~ l^w^. This models is a generalization of the multi- 
plicative process ( |To| ) and a practical implementation of a signal satisfying 
the multifractal representation ([)]). Studying the numerical process we 
found a "re-normalization" of the scaling exponents in space due the non- 
trivial interplay between multipliers re-generation and time evolution. We 
have shown that in the case of a passive scalar this effect is not present. 
Furthermore, we have clearly connected the asymptotic gap, observed on 
fusion rules, with the intermittent scaling exponents. This finding will be 
exploited in a forthcoming paper Q to build up a stochastic closure to 
compute anomalous scaling exponents for shell models of turbulence. 

We acknowledge G. Boffetta and A. Celani for many discussion in a 
early stage of this work. 
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Figure 1: Time behaviour of w n for n = 2, 4, 6 
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Figure 2: Log-log plot of measured scaling of Ss(n) vs. n for a time de- 
pendent binomial random multiplicative process. The slope fitted slope 
£ 3 ~ —1.36 is "renormalized" with respect to the "bare" value (3 = 1. 
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Figure 3: Behaviour of rcnormalizcd 0?(3) for log- normal distribution as a 
function of a 2 (see text). 
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Figure 4: Behaviour of Cr(p) f° r a — 0-03 (x) as compared to bare expo- 
nents Co(f) (+)■ Dashed line is the prediction (pif). 
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Figure 5: Time correlations at the same scale, (w3(t)ws(0)) , and at differ- 
ent scale, (u>9(t)u>3(0)), as a function of the time lag r. 
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Figure 6: The time scale B, computed by the expression 
(w n (t + T)w n {t)) ~ exp (— B|t|), for different value of n. The scal- 
ing constrain (ph should correspond to B ~ const as observed. 
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Figure 7: Log-log plot of (6^) vs n, for the passive scalar as compared to 
the theoretical prediction with slope —1. 
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Figure 9: Behaviour of Cg^t) as a function of time. 
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Figure 10: Behaviour of C 2 , 2 (oo) vs. ( R (4) - 2( R (2). 
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Figure 11: Behaviour of f(m) as a function of m for a log- normal distribu- 
tion and two different values of a, namely a = 0.03 (crosses) and a = 0.05 
(full line). 
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Figure 12: Plot ol the fusion rule coefficient G(m) , as defined in the text, 
for the passive scalar and for m = — 1 , . . 1 1 . 



